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In this paper we show that the homogeneity of the Pfaffian form representing the infinitesimal heat 
exchanged reversibly 8Q rev by a thermodynamic system allows to find immediately and explicitly 
an integrating factor. An interesting bridge between Caratheodory's approach to thermodynamics 
and Gibbs' approach is established. 
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One can distinguish in general between two main formal approaches to thermodynamics. On one side there is the 
approach due to Caratheodory, which is based on the integrability of the Pfaffian form 6Q re v [1—12] and represents the 
development of the line of thought which starts with Clausius and Kelvin. On the other side one finds the approach 
due to Gibbs, in which S is postulated to be an extensive concave function of the extensive variables [13-15]; further 
developments of this approach, from the point of view of its axiomatization, find their latest realization in Ref. [16]. 
Our aim in this paper is not to propose an axiomatics for thermodynamics from first principles, but simply to 
• discuss a formalism which represents a straightforward link between Caratheodory's approach and Gibbs' approach 
to thermodynamics and to propose a set of constructive assumptions which lead to standard thermodynamics. 
We know that Caratheodory's postulate of adiabatic inaccessibility ensures the integrability of the Pfaffian form SQ rev 
[3-7,9,10]; in particular the so-called metrical entropy S and the so-called absolute temperature T exist and 

(N ! 6Q rev = T dS. 

O ■ 

We choose to work in the framework of Caratheodory's approach, thus we assume the integrability of the Pfaffian 
form SQ rev . Moreover, we apply this formalism to (maybe multi-component) homogeneous simple systems, and choose 
the extensive variables (U, V,X, . . . ,X n ) as independent coordinates in the equilibrium thermodynamics space [for 
i-^h | "extensive" we mean "positively homogeneous functions of degree one"]. Notice that these variables characterize 
Gibbs' approach in the entropy representation, they characterize the so-called Gibbs space [15]; moreover they can be 
considered as the natural variables to be concerned with in the framework of an approach where the entropy and the 
• • ■ temperature have not yet been introduced. In the following, the main assumptions of our construction are indicated 
" as hl),h2),h3),h4),h5),h6),h7),h8),h9),hl0) in the text. We start from hi): 

x : 

hi) In our approach the homogeneity of the system is translated into the homogeneity of degree one of the integrable 
Pfaffian form 

5Q rev =dU +pdV £idX\ (1) 

i 

This means that, under the rescaling (U, V, X 1 , . . . , X n ) i— > (At/, XV, XX 1 , . . . , XX n ) one find 6Q rev ^ X SQ rev . See 
the following section for more details, p, £i, . . . , £ n are intensive variables. 

From a physical point of view, the homogeneity of SQ rev can be phenomenologically inferred by the observation 
of the extensivity of £7, V, X , . . . , X n and the intensivity of p, £x, • • • ; £ n m the case of standard thermodynamic sys- 
tems. Our main point is the following. Because of a symmetry of the Pfaffian form for homogeneous systems, we are 
able to find an explicit integrating factor for 5Q rev and to construct a potential which is identified with the entropy, 
in particular, with the entropic fundamental equation of Gibbs. In particular, both the entropy and the temperature 
appear as derived quantities. More general constructions for thermodynamics under the hypothesis of a non-trivial 
symmetry for SQ re v are studied in [17,18]. 
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II. HOMOGENEOUS PFAFFIAN FORMS 



In the following, we assume that the differential 1-form u) is at least C 1 in the domain V. We define il k (T>) as the 
set of all the k-forms (with a suitable degree of smoothness) defined on T>. We give below some relevant definitions. 
By "symmetry" for a differential 1-form ui we mean that, if Z is a vector field and Lz ■ ri 1 (P) — > SI 1 (2?) is the Lie 
derivative along it, we have 

L z wAw = 0. (2) 

The interested reader is referred to Rcf. [19] and in particular also to Rcf. [20], where the notions of homogeneous 
integrable Pfaffian forms and symmetry are found, and the corresponding integrating factor is deduced. 
The definition of homogeneous 1-form we give below is the same as the one appearing in [21]. An homogeneous 1-form 
of degree k + 1 

n 

uj(x\...,x n )=J2 0Ji{x\...,x n ) dx\ (3) 

i=l 

is, by definition, such that the coefficients uj i (x 1 , . . . , x n ) are homogeneous functions of degree k: 

0Ji(X x 1 , . . . , A x n ) = \ k oji(x\ . . . , x n ) V i = 1 . . . n. (4) 

In the case of an homogeneous 1-form the so-called radial vector field (which could be called also "Liouville field" 
[21]) 

Y = y (5) 

is a symmetry for u). In fact, if x stays for x 1 , . . . , x n , one has that the flow of the vector field Y is (f)t{x) = x exp(t), 
thus 

L Y lu = (k + l)uj. (6) 

If the form is integrable and co(Y) ^ 0, then oj(Y) is an integrating factor for ui. This follows from the integrability 
condition 

u A duj = (7) 

and from the definition of "symmetry". The condition u)(X) ^ for a generic symmetry generator X is equivalent to 
the requirement that the symmetry is nontrivial according to the definition of Ref. [19], i.e., the vector field X does 
not belong to the distribution P of codimension one associated with the kernel of w. This means that the symmetry 
"shuffles" the leaves of the foliation associated with the integrable Pfaffian form u) but it leaves the foliation itself 
invariant (instead a trivial symmetry preserves each leaf) [19]. From a geometrical point of view, a trivial symmetry 
is tangent to the foliation, a nontrivial one is transverse to the foliation. Note that in Rcf. [20] a homogeneous Pfaffian 
form such that ui(Y) ^ is called non-dicritical. We sketch here the proof that uj(Y) is an integrating factor for 
the sake of completeness; an equivalent but more elementary proof, based on a "by hand" verification, is found in 
Appendix A. / is an integrating factor for u> if 

-(7)-°- 

that is 

/ div - df A lu = 0. 

We verify that u)(Y) satisfies the latter equation: 

Lu(Y)du - (duj(Y)) Aoj = 
(iY(u))du) — ((d«y)(w)) Auj = 
(iy{i-o))duj — (LyiLi) Awl (iy{duj)) A ui = 
iy{oj A duj) = 

where we have introduced the standard contraction operator iy : SI 1 (2?) — > fl° (£>) and we have used (2), (7) and also 
standard identities (duj(Y) = (diy)(uj) = (Ly—iyd)(uj) = Ly ui — iy(duj) and iy (id Aduj) = (iy(Lo))dLo+((iyd)(uj))Abj). 
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III. THERMODYNAMIC ENTROPY REVISITED 



In the framework of thermodynamics of homogeneous systems, the domain V of 6Q rev is assumed to be an open 
convex set (in the discussion of the third law the boundary T = is introduced) . Convexity of T> is related to the 
request of a concave entropy (see the following subsections) . From a mathematical point of view, homogeneity forces 
the domain V to be closed under multiplication by a positive real scalar A, i.e. for each (U, V, A 1 , . . . , X n ) e P it has 
to hold (At/, XV, XX 1 , . . . , AA") G T>. This means that T> has to be a cone, thus it has to be a convex cone. Given a 
convex set C, a convex cone can be easily constructed, in fact the set 
K c = {{XU,XV,XX 1 ,...,XX n )\(U,V,X { ,...,X n ) eC,X > 0} 

is the smallest convex cone containing C and it has the property to be closed under addition (see [22], pp. 13-14). 
Notice that V is also a C°° differentiable manifold. 

The Pfaffian form SQ rev is assumed to be of class at least C 1 (X ) ), which is a rather general condition for the validity 
of Frobenius' theorem. [According to Rcf. [16], a more general setting should be allowed, where the intensive variables 
are only locally Lipschitz-continuous. This requirement, which is physically well-grounded [16], would require a more 
general setting for Frobenius' theorem. In fact, the Pfaffian form 5Q rev would be locally Lipschitz-continuous, which 
implies that d(8Q rev ) exists only almost everywhere, thus also the integrability condition is defined only a.e. We 
don't deal with this problem herein]. Then the intensive variables p(U, V, X 1 , . . . , X n ) and S,i(U, V, X 1 , . . . , X n ), for 
i = 1, ... ,n, belong to C 1 (I?) (at least); they are supposed to be known functions of the extensive variables. 
The Liouville operator is 



h2) We require that the homogeneity symmetry of SQ rev is nontrivial. 



Then, an integrating factor for 5Q rev can be immediately deduced from the homogeneity property, and it is given by 

f(U,V,X\...,X n )=i Y (u) = U+pV 6*^0. (9) 



As a consequence, 

SQ r 



f 

is a closed and, moreover, exact 1-form (being the domain simply connected). Thus 



(10) 



where S is a potential for the form. The 1-form dS is an homogeneous function of degree zero, as it follows from (11); 
in fact, the integrating factor / is an homogeneous function of degree one: 

Y f = Y iy(u>) = lyLyiiJ = iyUJ = f 

(the Cartan formula Lxiy — Lyix = i[x,Y] holds), and the infinitesimal heat exchanged reversibly is an homogeneous 
function of degree one. The exact Pfaffian form in (11) has coefficients (1//, p/f, — • • • , — £«//) which arc 
homogeneous functions of degree —1. Then S can be found only by quadratures (cf. [23], p. 16-20). One has 

S(U, V, X 1 , . . . , X n ) - S(Uo, Vb, Xq, . . . , Aq 1 ) = jf ^^T-e/Ia^ (12) 

where (C/ , Vo, Aq , . . . , Xq) is a reference state and T is any reversible path connecting the reference state to the state 
(U, V, A 1 , . . . , A") of interest [a path is a oriented curve in the thermodynamic domain which is at least piecewise C 1 ; 
one could also require, without loss of generality that the path is also simple, that is, non self-intersecting. A path, 
being defined along equilibrium states of the thermodynamic manifold, is, as a consequence, reversible in the sense of 
Refs. [7,14]]. 
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h3) We require that the thermodynamic foliation is defined by the leaves S = const, everywhere in D. 

As a consequence, S is required to be a well-defined function for each state belonging to the thermodynamic do- 
main T>. This assumption is natural in our framework; its necessity emerges in light of the discussion of the following 
subsections, in particular of subsections IIIC and HID. 

h4) We require that the integrating factor / is non-negative 

The non-triviality condition / ^ is enforced by requiring that / is of definite sign [one could also ask for / < 0, the 
choice / > is made on a conventional basis] . The relation of hypothesis h4) with thermodynamics is evident from 
the following subsection, being related to the non-negative definiteness of the absolute temperature. 



A. extensive entropy 

The potential S is related to the logarithm of a positive definite extensive potential H 

S-S = log (J^j , (13) 

where Y H = H and H = H (Uq, Vq, X„ , . . . , ) is the value of H at the reference state (Uo, Vq, Xq , . . . , Xft). 
Analogously, S = S(U , V , A \ Xff). 
In fact, one has 

dL Y S = L Y dS = L Y j =0, (14) 

being L Y (uj/f) = (L Y u))/f - (lo If) (Y f) = 0. Then 

LyS = YS = q (15) 
where q is a constant. We can, without loss of generality, define another positive definite function H such that 

S = \og(H). (16) 
Then we get that (15) is equivalent to the following equation for H: 

YH = qH, (17) 
which means that H is an homogeneous function of degree q. One has 

- dH lo , . 

dS = — = -; (18) 

as a consequence, 

dH=juj (19) 

which is an homogeneous of degree q closed Pfaffian form. Let us assume that q ^ 0, i.e., that H is not intensive. 
Then, from (A21) in Appendix A it follows 

H = 1 j i Y (u) = ^ H, (20) 

which implies q = 1. The case q — is treated in Appendix B, and it is excluded, H cannot be intensive. See also 
[17], where a more general proof is given. 

The leaves S — const, of the thermodynamic foliation coincide then with the leaves H — const, and H appears to 
play a privileged role in the construction of the thermodynamic foliation. 
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B. metrical entropy S 



We know that there exists a function H such that 

m v,#,..., xn - s { uo, Vo ,xi ...,xs)= io g < 2i > 

where H is an extensive function. We have 

SQrev =fdS =L dH, (22) 
H 

where f/H appears as a new integrating factor and H as extensive potential. By symmetry, the function H is 

suggested to be the thermodynamic potential one is looking for, i.e. the so-called metrical entropy [3], because one 
can write 

jj dH = dU + P dV -J2 & dxi ( 23 ) 

where f/H is intensive as each coefficient of the Pfaffian form is, and H is extensive as each independent thermo- 
dynamic variable U, V, X 1 , . . . , X n is. This identification of H with the metrical entropy is correct. We know, from 
the standard approach [3-7], that a metrical entropy has the property to be additive. In particular, the entropy S 
of a system composed by two subsystems in thermal equilibrium is S — Si + S*2. Moreover, it is expected that the 
entropy of a system composed by n replicas of a subsystem E is S — n So, that is, it is extensive. We just have 
the extensivity property for H. Moreover, also on a statistical mechanical ground, S is required to be non-negative. 
Then, we can introduce on this ground the following assumption: 



h5) We require that the metrical entropy S is extensive. 



We can easily show that H is a metrical entropy because, up to a multiplicative constant, there exists only one extensive 
potential S and only one intensive integrating factor T such that 5Q rev = T dS; as a consequence, f/H = (dH / dU)^ 1 
plays the role of absolute temperature. The proof is the following. We have obtained 

uj=g dH, (24) 

where g = f/H is homogeneous of degree zero. One may wonder if g and H arc unique. By introducing another 
function G(H) we can obtain 

"stm" 3 - (25) 

We want again the new integrating factor g = g/(dG/dH) to be homogeneous of degree zero and G to be homogeneous 
of degree one. Then we have 

Y G = G (26) 

and 



Then 



which means 



C-^H (28, 



G = 7 H (29) 

where 7 is a constant. As a consequence, one finds 
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1 8G dH 1 

The entropy H and the absolute temperature g are then unique apart from a scale factor 7 which can be fixed by 
fixing the absolute temperature scale [3] . S plays instead the role of empirical entropy [3] . 
We have found that, apart from a scale factor, 

H = S, (31) 
d — = - - (32) 
The latter equality, together with the standard choice for T to be non- negative, justifies h4). 



C. corroboration from Gibbs' approach 

S can also be easily related to the standard definition of entropy also by simply appealing to Gibbs' approach. 
In fact, it is sufficient to consider the fundamental equation of thermodynamics in the entropy representation for an 
homogeneous system [14], which implies 

TS = U+pV-^^X i ; (33) 

i 

then 

f = TS (34) 

this shows that / cannot identically vanish in thermodynamics, thus, a posteriori, the integrablc Pfaffian form SQ rev 
has to be non-dicritical, and the homogeneity symmetry has to be non- trivial (or transversal). Moreover, 

dS=^-=dlog(S). (35) 

Of course, both (34) and (35) could be easily found at first sight from (9), but we have shown that on the ground of 
simple assumptions one can recover S and T as derived quantities without referring a priori to Gibbs' fundamental 
equation in the entropy representation. 
We find that 

^ ~ = f -• (36) 
So & Jr f 

As a consequence, 
and 

S = S exp(S - S ) (38) 
= So exp(^ j). (39) 

Notice that, a priori, one should write dS = dlog(\S\), and insert an absolute value in the argument of the logarithm 
in (37). In fact, from the point of view of classical thermodynamics, there is no constraint for S to be non-negative. 
Nevertheless, the left side of (36) has to be of definite sign, contrarily S would not be well-defined where S = 
(the left side of (36) could still be defined as a principal value integral) and the thermodynamic foliation defined by 
S would not be in a one-to-one correspondence with the one defined by S, because the map S ^ S would not be 
injective. 

From (37) we can also deduce that S is not an homogeneous function, in fact a differentiable homogeneous func- 
tion g(U, V, X 1 , . . . , X n ) of degree a should satisfy Y g — a g. For S we get 

L Y S = Y S = 1. (40) 
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D. positivity and concavity of entropy 



It is interesting to notice that the thermodynamic entropy S is positive. There is no need of restricting the range 
of the metrical entropy S to positive values, because the thermodynamic potential H, which is realized to be a 
metrical entropy, is, by definition, positive definite. Of course, this agrees with the statistical mechanical approach. 
The classical ideal gas violates this property, in fact its entropy becomes negative and goes to — oo as T — > 0; but 
this behavior being corrected at low temperatures by quantum mechanics. [The positivity of S can be obtained in 
the framework we sketched in the previous subsection by recalling that, under the assumption h3), S is required to 
be well-defined on each thermodynamic state belonging to the thermodynamic domain where the thermodynamic 
foliation is constructed; see the discussion in the previous subsection]. States such that S = are singular and can 
be allowed to belong at most to the boundary of the thermodynamic domain, as we show in the following subsection. 

As far as the concavity of the entropy is concerned, we limit ourselves to discuss how this fundamental stability 
requirement which is introduced in Gibbs' approach constraints the Pfaffian form SQ rev . Concavity of S ensures that 
the system is thermodynamically stable [14]. We introduce then the assumption 



h6) We require that the metrical entropy S is concave. 



Notice that the requirement of superadditivity for S would be equivalent, in fact superadditivity and homogene- 
ity of S imply concavity [24,25]. Cf. also [17]. 

First, we underline that the concavity of S is not equivalent to the concavity of S. From one hand, the concavity of 
S is a sufficient condition for the concavity of S; in fact, the logarithm of a concave function is a concave function (if 
g is a concave function and is a concave non-decreasing function, then <j> o g is a concave function. The logarithm 
is a concave non-decreasing function) . The concavity of S is not a necessary condition for the concavity of S (if one 



considers S = exp(— x ), which is not concave on IR but only for 2 x — 1 < 0, one still obtains S 



which is 



concave on IR). On the other hand, the concavity of S does not ensure the concavity of S, because, in general is not 
true that the exponential of a concave function is concave (e.g., let us consider again S = —x 2 , which is a concave 
function on IR. The function S = exp(— x 2 ) is not concave on the whole real line). 

The concavity of S can be obtained by imposing suitable conditions on the coefficients of the Pfaffian form 5Q rev . 
These conditions can be deduced from the standard ones on S [14]. Let us consider for simplicity S(U, V, N). We get 



dS_ 
dU 
dS_ 
dV 
dS_ 
dU 



= S 



= S 



The Hessian matrix for S is easily deduced to be 



D 2 S = 



S 



1 dU 



P dV 



r dN 



1 

7 
p 
f 

f ' 



P dV 



(41) 
(42) 
(43) 



r)V f av 



' ON 
-VP + f^k-P§h 



-HP 



f ®1L — n d/ 
■I dN y dN 



J dN " dN 



Concavity requires that all principal minors of odd order should be negative and all principal minors of even order 
should be positive. The Hessian determinant is zero because of the homogeneity of S. Then, in order to obtain 
a concave S one has simply to impose on the homogeneous integrable Pfaffian SQ rev the following (necessary and 
sufficient) conditions involving the coefficients of the Pfaffian form and the integrating factor /: 



df 

1 -dU <0 

n d f u 2 , r dp 9f 

{1 -du )ip +f W- p dv ) 



( 9 A 



> 0. 



(44) 
(45) 



Notice that 1 — df/dU < amounts simply to dT/dU > 0, which is simply the positivity of the heat capacity at 
constant V, N. 
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E. conditions for S and / to be globally denned 



The homogeneity of the integrable non-dicritical Pfafhan form 6Q rev allows to find an integrating factor / whose ex- 
pression holds globally by construction (by hypothesis, p(U, V, X 1 , X n ) G C 1 ^) and &(t/, V, X 1 , . . . , X n ) G C^D) 
for alH = 1, . . . , n). A priori it is not possible to ensure that also S is defined everywhere, because singularities can arise 
where / = 0, i.e., if the set of zeroes of the integrating factor Z(f) = {([/, V, X\ ... , X n )\ f(U, V,X\..., X n ) = 0} 
is non-empty. In fact, where / vanishes S can be not well-defined, as it is clarified in the following, in contrast with 
our requirement that the thermodynamic foliation is described by the potential S in the convex set V. 

1. the set Z(f) 

We know that f = TS> 0, and obviously Z(f) = Z(T) U Z(S), where Z(T), Z(S) arc the sets where T, S vanish 
respectively. Surely S is not well-defined for each state belonging to Z(S), because of (37). Thus, one has to require 
that Z(S) belongs at most to the boundary of V. This requirement is corroborated by the analysis of convex functions 
[26], in fact / = — S is a convex function which attains its maximum value in Z(S) if Z(S) ^ (recall that S > 
in our framework), and, in order that / is non-constant, it is necessary that Z(S) is contained in the boundary of 
the convex domain V [cf. also [27], thm. C p. 124]. Z(T) is contained in the boundary as well (as it is evident if T 
appears as independent variable). Then also Z(f) is contained in the boundary of the thermodynamic domain. 
Moreover, 

h7) We require that Z(S) C Z(T). 

In the following discussion, we indicate with z collectively the independent thermodynamic variables. The assump- 
tion h7) is introduced because any state z such that T z > and S(z) — should belong to the boundary of the 
thermodynamic domain and should have the peculiar property to allow the system only to absorb heat along any 
non-adiabatic path j z starting from z in a neighborhood W z of z. In fact, let us define the heat capacity along a path 
7 : [T , T] — > T> which does not include isothermal sub-paths: 

C 7 (T)^T (g) = W (7). (46) 

If 7 Z is a path starting from the state z, then 

S (y) = J T -jrC 7t (T) (47) 

should be positive for any state y non isoentropic to z in W z , because S(y) > S(z) = 0, which is possible only for heat 
absorption (in fact, C 7z (T) < would be allowed for states such that T z <T yi which would imply heat absorption, 
and C lz (T) > would be allowed for states such that T z > T y ). [Heat absorption should occur also for any isothermal 
path starting from the state z, in fact Q — T z AS along an isothermal path and AS > 0]. Thermal contact with a 
colder body at T < T z should allow an heat flow outgoing from the system (see also the discussion in Ref. [26] ) . Then, 
no quasi-static approximation of such a thermal contact can be allowed, no matter how near to T z the temperature of 
the colder body could be, because, at least in a neighborhood of T z , the system could only absorb heat. This behavior 
can be considered pathological, and the occurrence of the absolute minimum S — of the thermodynamic entropy 
at T > is refused in the framework of standard thermodynamics. As a consequence of the rejection of states with 
S = at T > one finds Z(f) = Z(T). Notice that in the classical ideal gas case, where S is allowed to become 
negative, one finds that Z(f) D Z(T) because T vanishes for U = and S vanishes before the hypersurface U = is 
reached. 

Notice that h7) is automatically implemented if the boundary T = is described by a (maybe even smooth) 
function U = biV^X 1 , . . . , X n ) and conditions ensuring the continuity of S at T = are allowed [26]. The function 
b(V, X 1 , . . . , X n ) can be construed as a ground-state energy. h7) can be substituted by the stronger assumption 
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h7) We require that T = is described by U = X 1 , . . . , X n ), where b is an extensive convex function defined on a 
convex cone JCb, and that the domain T> U dV coincides with the epigraph of b: 

Z> U <9£> = cpi(6) = {(U,V,X\...,X n ) | (V,X\...,X n ) eJC b ,U> b{V, X 1 , . . . , X n )}. (48) 

As a consequence, by defining the extensive coordinate B = U — b(V, X 1 , . . . , X n ) > it is easy to show that 
Z(S) C Z(T), because for an everywhere continuous entropy it holds 

S(B,V,X 1 ,...,X n ) = S(0,V,X 1 ,...,X n ) + J* _ (49) 

where S(0, V, X 1 , . . . , X n ) is the value attained by S at B = by continuity [all the matematical properties ensuring 
the existence of the improper integral dY l/T are implicitly assumed]. See [26] for further details. Then we have 
Z(f) = Z{T). Notice that, if b = 0, then V U dV = K. x IR^, where K, 3 (V, X 1 , . . . , X n ) is a convex cone. 

Then S and S are defined everywhere, because SQ rev /f is a closed Pfaffian form defined everywhere on V. This 
solves the problem of ensuring the global existence of the integrating factor and of the potential, in particular, the 
entropy and the temperature are defined globally (about this problem in the frame of Caratheodory approach, cf. 
[9,10])! 

We further introduce the following assumption: 
h8) We require that to each level set S = const, corresponds a unique leaf. 

In general, S is a submersion whose leaves are the connected components of 5 _1 (c), where c G IR + is a constant. 
From a physical point of view, it is expected that each isoentropic surface S = c = const, is path-connected, in fact, 
given a state X and a state Y lying on the same isoentropic surface of an homogeneous thermodynamic system, it is 
physically mandatory that there exists an adiabatic reversible transformation [i.e. an adiabatic path] which connects 
X and Y. Thus, assumption h8) has to be ensured. In Appendix C we show a possible implementation of assumption 
h8). 

It can be noted that the conditions lo(Y) ^ and lu(Y) = only on the boundary of the thermodynamic domain 
ensure that, given a fixed reference point zq = (U~o, Vo,Xq, . . . ,Xq), the integral curve 7z (£) of the vector field Y 
meets each leaf S = const, of the foliation if h8) is implemented, in fact S(T>) = (0, oo) and it holds 

( 7 ; (t)S)(U , V , Xl ...,*?) = S(exp(i) U , exp(i) V , exp(t) A \ . . . , X ") = exp(i) S(U , V a ,Xl . . . , X%). (50) 

Moreover, there exists only one intersection with each leaf. The given curve is then transverse with respect to the 
foliation. 

Another physically well-grounded assumption is 
h9) We require that J|r is positive. 

This requirement means simply that the pressure p is positive definite. 

F. the problem of T = 

This kind of approach is very interesting also from the point of view of the third law of thermodynamics [3-5,7,28]. 
A more complete analysis of this topic is the subject of Ref. [26,29]. We limit ourselves to a short summary of some 
points which have an evident link with the subject of this paper, referring the reader to Ref. [26,29] for details, for a 
general discussion and for also for complete references to the literature on the third law. 

We assume that the T = is a connected hypersurface (see also [26] ) , which coincides with the adiabatic boundary 
SQrev — of the thermodynamic domain, i.e. 
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hlO) We require that T = is a connected integral manifold of 5Q rev . 

This means that, if h7) is implemented, U = b(V, X 1 , . . . , X n ) has to be a solution of the Mayer-Lie system which is 
equivalent to the Pfaffian equation w = (see [12]), i.e. the following conditions 

^L = - p (U,V,X\...,X n ) = § F (51) 

8U db 

— = ^U,V,X\...,X n ) = — fori = l,...,n (52) 
have to be implemented. Then the function b has to be at least of class C 1 (/Cb). 

Requirement hlO) is not necessary if h7) is implemented and / and to are C 1 everywhere (then b has to be of class 
C 2 (/Cb)). In fact, from / dio = df Au> one finds that the pull-back of the map Fb : V — > dV, where F b *(f) = f o Ft = 0, 
implements = F b *(df A w) = F b *(df) A F b *(co). This implies that there exist a real function /i ^ defined in 92? such 
that F b *(ui) = h F b *(df) and F b (df) = because F b *(f) = 0. Cf. also [20]. Under more general conditions on / and u 
it can be non-trivial that / = is an integral manifold. 

The use T as explicit independent coordinate is not suitable for studying the nature of the integral submanifold 
T = 0, in fact the Pfaffian form 5Q rev seems to be singular in T = (the point is that the map U ^ T is not a 
diffcomorphism in T = 0) . In order to avoid this problem the explicit use of U or B [or of another regular coordinate] 
is mandatory. We then use B explicitly See also the discussion in [26,29]. Outside T — [i.e. B — 0] there exists 
a thermodynamic foliation whose leaves are the hypersurfaces 5 = const. The nature of T — is ambiguous, in the 
sense that it can correspond to a leaf of the thermodynamic foliation, the leaf T = 0, or it simply can correspond to 
an integral manifold of SQ rev which is intersected by other integral manifolds of SQ rev . In both cases, S : V — ► IR+ 
represents an at least C 2 (V) submersion which generates the thermodynamic foliation at T > 0. We can show that 
the validity of the third law of thermodynamics can be related with the geometric nature of the integral submanifold 
T = 0. 

In fact, under the hypothesis that the entropy 5 is concave and continuous at the surface T = 0, in the sense that 
there exists the limit as B — > + of S for any finite value of the parameters V, X 1 , . . . , X n [notice that the theory of 
convex functions allows to define S at T = 0, i.e. when U = b [26]], it can be shown that T = corresponds to a 
special leaf of the thermodynamic foliation generated by the Pfaffian form SQ rev if, and only if, the entropic version 
of the third law 

lim 5 = (53) 

holds. We notice that Planck's restatement of the third law (53) is trivially mandatory for an homogeneous 
system, in fact lim B ^ + S(B, V, X 1 , . . . , X n ) = So, with 5 > a positive constant, implies also that So = 
lim B ^ 0+ S(XB, XV, XX 1 ,..., \X n ) = X lim B ^ + S(B, V, X 1 , . . . , X n ) = XS , which is possible only for S = 0. It 
can also be noticed that the positivity and the concavity of S forces S to be finite as B — > [26]. In fact, the behavior 
of the convex function I = — S at the boundary has to be such that 

liminf I{x) > -oo (54) 

for any x = (U, V, X 1 , . . . , X n ) which converges to xo = (Uq, Vq, Xq, . . . , Xff) belonging to the boundary of the convex 
domain [cf. problem F p. 95 of Ref. [27]]. Then a positive S has to be finite as T — > 0. 

If (53) holds, then Z(S) = Z(T) and no leaf S = const, can intersect T = 0, in fact only the isoentropic surface 
5 = could intersect T = 0. If (53) is violated, and lim B ^ 0+ S(B, V, X 1 , . . . , X n ) = S(0,V, X 1 , . . . , X n ), then 
Z(S) C Z(T) and T = is not a leaf; in fact, one can uniquely define the entropy to be 5(0, V, X 1 , . . . , X n ) at 
(B = 0, V, X 1 , . . . , X n ), and a peculiar intersection of adiabatic surfaces occurs [26]. The main point is that, if (53) is 
violated, a state belonging to the set Z(T) can also belong to an isoentropic surface intersecting T = 0, and there are 
two adiabatic paths starting from that state, one contained in T = and the other contained in the aforementioned 
isoentropic surface. This happens as a consequence of the fact that the surface 5 = S(0,V, X 1 , . . . , X n ) = const. 
> 0, if (0, V, X 1 , . . . , X n ) is not a local minimum of 5, necessarily intersects the graph of 5 along a codimension one 
submanifold which is connected to the submanifold T = 0. [The proof is found in [26]; see also [29]]. This means 
that the integral manifold 5 = 5(0, V, X 1 , . . . , X n ) = const. > reaches T = and is actually tangent to it [26]. 
Notice that, given such a point (0, V, X 1 , . . . , X n ), then each point (0, XV, X X 1 , . . . , X X n ) belonging to the cone of 
(0, V, X 1 , . . . , X n ) is such that 5(0, A V, X X 1 , . . . , X X n ) = X 5(0, V, X 1 , . . . , X n ) and it is not a local minimum as 
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well [the proof is trivial]. Then, each (would-be) leaf S = const, intersects T — 0. One does not obtains a well-defined 
thermodynamic foliation if T — is included, actually a foliation can be obtained in this case only if zero-measure 
set represented by the surface T = is excised from the thermodynamic domain. Contrarily, one obtains a foliation 
except for a zero-measure set represented by the surface T = 0, i.e. the inner part of the thermodynamic domain 
is foliated into surfaces S = const., and locally the adiabatic inaccessibility holds, but these (would-be) leaves are 
actually connected by the adiabatic surface T = and the adiabatic inaccessibility is violated, even if only along 
special curves reaching T = [26]. Note that the failure of the Lipschitz property for ordinary differential equation 
for adiabatic curves starting at T = is related to the singular behavior of the thermodynamic foliation at T = 
when (53) is violated [26]. Note also that, on the ground of assumption h8), when (53) is violated one has to impose 
that any point (0, V, X 1 , ... , X n ) which is a local minimum has to be also a global minimum (contrarily, h8) would 
be violated). 

If, instead, T = is a leaf, then no leaf S = const, is allowed to reach T = 0. Leaves S — const, can at most 
asymptotically approach T = 0, but no one can be extended therein. (53) has to hold for a continuous S because the 
limit of S as T — > cannot depend on the parameters V, X 1 , . . . , X n (otherwise T = is not a leaf) [26,29]. 

Note that also the approach discussed herein is able to make immediately evident the necessity of the discussion 
of the states which belong to Z(f ), i.e., as a consequence of h7), of the states which belong to Z(T). The third law 
has to be discussed in the framework of thermodynamic formalism, and the validity of (53) emerges as a "regularity" 
condition on the boundary T = for the differential equation SQ rev — of thermodynamics, in particular it appears 
as a condition which has to be imposed if a foliation of the whole thermodynamic domain (including T = 0) has to 
be obtained. 



In the framework of our approach, (53) holds if and only if, whichever path 7 is chosen in approaching T = and 
having a point belonging to T = as its final point, the integral of the form 5Q rev / / diverges to —00. The condition 
is clearly necessary, because (53) implies the divergence of (37); its sufficiency is also evident, being equivalent to 
log(S') — > —00. Notice that this condition on the integral of ui/f docs not require to define the entropy at T = 0. 

Moreover, it is easy to show that a sufficient condition for (N) to hold is u> € C 1 everywhere [i.e. also on the 
boundary T = 0], which implies / G C 1 everywhere. In fact, one has 

^L-l + S^ (55) 

which has to be finite at T — 0; it holds dT/dU = l/C x ^.... t x n + 1 , where Cxi,...,x™+ 1 is the standard heat capacity, 
which has to vanish in the limit as T — > because S has to be finite in that limit. As a consequence, dT/dU = 
l/C x \...,x-+i -> 00 as T 0+, thus (55) is finite only if S as T 0+. 
Further details and conditions are the subject of forthcoming papers [26,29]. 



G. a note on the reference state 



The reference state (Uo, Vo, Xq , . . . , Xff ) has to satisfy a constraint which is trivial but its discussion can be physically 
interesting. From the point of view of thermodynamics, given a system described by three variables (U, V, N), reference 
states like (0,0,0) or (0, V, 0) are meaningless. In fact, they correspond to the absence of the system under study. 
Statistical mechanics provides a criterion for obtaining a meaningful thermodynamic description: there should be a 
statistically relevant number of elementary constituents (atoms, molecules, particles) of the system. In particular, the 
thermodynamic limit is well-known to be the tool allowing to find out thermodynamics from statistical mechanics. 
[See also [17]]. N — correspond to the absence of constituents. Without necessarily referring to the above statistical 
mechanical suggestion, one may wonder when a thermodynamic description starts being available. The system 
should be such that thermodynamic fluctuations allow to define meaningful macroscopic variables; moreover, the 
system under study should be realized to be involved with appropriate thermodynamic properties. In other terms, 
a "preparation" procedure of the system, in such a way that a macroscopic description becomes available and a 
set of relevant thermodynamic variables can be identified, should be considered as a preliminary condition to the 
thermodynamic description. 
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H. fundamental equation and homogeneous PfafHan form 



It is known that the fundamental equation contains all the relevant thermodynamic information about the system 
by construction [14]. In order to determine the fundamental equation S by integration of the Pfaffian form 5Q revl 
it is necessary to consider the Pfaffian form SQ rev representing a sufficiently general thermodynamic transformation 
which can be implemented with the system under study. Before clarifying what we mean by "sufficiently general" , let 
us discuss some examples. In the case of the classical ideal gas, one should consider the system as open and let the 
extensive variables U, V, N to appear in SQ rev . As it is easy to show, by considering SQ rev with N = const, (closed 
system) leads to a wrong fundamental equation, to be compared with the correct one displayed in Ref. [14]. If a 
photon gas is studied, only two variables are necessary, U, V, whereas TV is not a good parameter [which is related to 
the fact that the number of particles is not a good observable for a massless particles gas] . These examples show that 
neglecting some variable one would find incorrect results. Moreover, the relevant extensive variables should be settled 
out by considering, for each system, its own macroscopic features. Let us suppose that a system can be consistently 
described, when e.g. insulated, by means of three extensive variables U, V, N appearing in the fundamental equation 
S(U, V, N). This does not exclude a priori the possibility to consider the system in a more general situation in which 
the system is under the influence e.g. of some external field which could give rise to new thermodynamic features 
[as a paramagnetic gas under the influence of an external magnetic field, one could introduce an external magnetic 
field B e and the magnetization M as in Ref. [14]]. Then one can allow an extension of the thermodynamic space, by 
enlarging the set of the extensive variables, e.g. by passing from U, V, N to U, V, N, X, where X is the new extensive 
variable which can assume the value X — and whose conjugate intensive variable is £. A new fundamental equation 
S(U, V, N, X) is obtained, which, by continuity, reduces to S(U, V, N) when X = 0: S(U, V,N,X = 0) = S{U, V, N). 
Of course, an enlarged set of integrability conditions should be verified [cf. eqn. (59) in the next section]. A consistency 
condition for such an extension to hold consists in requiring that £(J7, V, N, X = 0) = 0. This can be easily realized if 
one refer to the Massicu potentials which are obtained by performing a Legendre transformation of the fundamental 
equation [14]. If <j){\/T,p/T 1 /i/T, £/T) is the Massieu potential associated with S in which only intensive variables 
appear, then, in order to obtain a consistent reduction of thermodynamic degrees of freedom by posing X = 0, it is 
necessary to put £ = 0. 

Then "sufficiently general" means that in the expression of the Pfaffian form SQ rev one takes into account all the 
extensive variables which cannot be consistently set equal to zero. This kind of approach to the problem of the 
construction of the fundamental equation appears to be more realistic than the one requiring to refer to the Pfaffian 
form SQ rev corresponding to the most general transformation the system could be allowed to perform [the most 
general transformation is hardly known]. An explicit construction of this kind occurs in black hole thermodynamics, 
and it is discussed in Ref. [17]. Because of the presence of gravity, thermodynamics is no more homogeneous, but a 
quasi-homogeneous symmetry of the Pfaffian form SQ rev can still be identified [17]. We point out also that quasi- 
homogeneity is the natural symmetry for standard thermodynamics when one or more independent variables are 
intensive [30]. 



In the standard approach to homogeneous systems, in the so-called entropy representation, it is known that the 
knowledge of the fundamental equation is equivalent to the knowledge of the so-called state equations [14]. In 
particular, we have the following general scheme for reconstructing the fundamental equation of thermodynamics 
from the state equations [14]. This method is equivalent to the one we propose, where we know the intensive variables 
p, £i, . . . , £„ and we reconstruct S. Given a system described by n + 2 extensive variables, it is sufficient to know n + 1 
state equations, i.e. n+1 intensive functions (pressure, chemical potential,...) in order to reconstruct the (n + 2)-esime 
state equation within an additive constant, and, then, to reconstruct the fundamental equation. In fact, being S an 
homogeneous function of degree one of the extensive variables, one has 



moreover, by differentiating the above formula for S, the Gibbs-Duhem equation follows, which allows to reconstruct 
the (n + 2)-esime state equation from the other ones: 



IV. COMPARISON WITH THE STANDARD APPROACH 




(56) 




(57) 
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Here, we assume to know p, £1, . . . , £ n and we can recover T (within an integration constant) from the knowledge 
of p, £1, . . . , £„. Then we can reconstruct S. Of course, in general such a reconstruction of the (n + 2)-esime state 
equation is possible only under suitable integrability conditions to be satisfied by p, £1, . . . , £ n . These integrability 
conditions coincide with the integrability conditions of SQ rev . A Pfaffian form 

uj = V" LUi dx % (58) 

i 

is integrable if it satisfies the integrability condition w A dui = 0, that is, if, for any i, j, k 

( dujj duj k \ ( duj k dojA ( duji dwj\ 

hjk=Wi {d^-d^) +Wj W-d^) +LUk {d^-d^) =0 - (59) 

[See e.g. Ref. [3]]. These conditions have to be satisfied by (1). E.g., in the case of three variables U, V, N, one has 

dp d/j, dp d[i 

luvN ^m + dv + ^du- p du =0 - (60) 

These conditions ensure that SQ rev / f is an exact differential. Let us consider the Gibbs-Duhem equation in which T 
is unknown 

d log IrJ = - u + pv-Xt m t = 1 5 (61) 

the integrability conditions ensure that (61) has a right member which is an exact differential, as can be easily verified. 
In fact, let us consider 

SQrev _ dU+pdV-J2 t ZidX* _df _ Vdp-J2 t 

~r _ 7 "7 i ; (6) 

it is evident that dS = SQrev/ f holds if and only if 

Vdp-Y.i X 1 dUi 



f 

is the differential of a function. Moreover, from / = T S and from 



(63) 



dS= d 4- Vdp ~^ Xi ^ = d A (64) 
./ ./ 6 

the Gibbs-Duhem equation (61) follows. Note that the Gibbs-Duhem equation (61) allows to find log(l/T) within 
an additive constant, which means that T is determined within a multiplicative constant, and also S is determined 
within a multiplicative constant. In in this respect, the integration of uj/ f and the integration of (61) are equivalent. 
As a counter-example where the integrability conditions fail, let us consider SQ re v = du + pdV — fi dN, with p = U /V 
and \x = UV/N 2 . The integrability condition (60) is not satisfied, in fact Iuvn = U/N 2 ^ 0. One has 

/ 1 \ V - N V + N 2V 

dl ° g bJ = U(2N -V) dU+ 2N^V dV N(2N V) ^ 
= au dU + ay dV + dN. 

The right member is not an exact differential, because e.g. dyau = N/U(2N — V) 2 ^ djjav = 0. About the 
Gibbs-Duhem equations see also [17] and [30]. 

A. densities 

We introduce the so called densities. One defines density the ratio of the extensive independent variables with 
respect to one of them. For example, it is possible to define: 

u = U/V 
x l = X l /V. 
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The n + 1 densities so obtained are the physical degrees of freedom in the case of an homogeneous substance. Then 
the fundamental equation in the entropy representation becomes 

S(U,V,X\...,X n ) = V s(u,x\...,x n ), (65) 

where 

8{u,x\...,x n ) = ^u+^-^x i , (66) 

i 

where the intensive variables are now functions of u 7 x\ . . . , x n . The Gibbs-Duhem equation becomes 

ud (^) +rf (f)-£^(i) =0 ( 6? ) 
and now only intensive variables appear. 

In the approach by means of the homogeneous differential form SQ rev , the use of the densities leads easily to the 
following result: 

SQ rev = (u+p-J2 &x l )dV + V (du-J2 6 (68) 

i i 

Then, we get 

SQrev , dV 

-J— = + — , (69) 

where 

_ du - J2< 6 dx l ds 
w o = — ; : — j — — • (,'U) 

Notice that in the latter approach both s and T are unknown. All the intensive variables p, £i, . . . ,£„ are to be known 
if the entropy has to be determined, thus p has to be known in order that the fundamental equation can be recovered. 
This is completely in agreement with what happens in the standard approach, where only one intensive variable is 
allowed to be unknown [it can be reconstructed by means of the Gibbs-Duhem equation]. One may wonder why in 
the case of an homogeneous one-component system one can sometimes recover the fundamental equation S(U, V, N) 
without any knowledge of the chemical potential /x. The point is that, by passing to the densities with respect to N 
[indicated again with u, v], one has S — N s(u, v) and 

dU + pdV-iidN 1 , . 1 , , . , ,„ 

dS = £-= = r N (du + p(u, v) dv) + — [u + p(u, v)v- fj,(u, v)) dN (71) 

1 1 (u, V) 1 (u, V) 

= du+p(u,v) dv u + p(u,v)v - n(u,v) 

T(u,v) T(u,v) [ ' 

= N ds + s dN. (73) 

Then, it is possible to find s(u, v) by considering the system as closed if the temperature T(u, v) is known. The point 
is that it is legitimate to ignore fj, only if the dependence of the temperature on u, v is known. A typical example of 
such a system where T is known is represented by the classical ideal gas [14]. 



V. CONCLUSIONS 



In this paper we have shown that, by choosing the extensive variables (U, V, A 1 , . . . , A") as independent variables, 
the intcgrability and the homogeneity of 5Q rev allow to find explicitly and immediately an integrating factor for 
8Q rev . The relation between the potential S and the entropy S is straightforward. It is to be noted that the 
entropy S arising from this calculation corresponds to the fundamental equation in the entropy representation. Thus 
our approach allows to find a direct link between Caratheodory approach, based on differential forms, and Gibbs 
one, based on the postulate of the concavity and homogeneity of the entropy S(U, V, A 1 , . . . , A"). The role of the 
homogeneity as a symmetry for the system and as a tool for constructing the thermodynamic formalism has been 
enhanced. 
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APPENDIX A: INTEGRABLE HOMOGENEOUS PFAFFIAN FORM 



Here we sketch a more elementary proof for some statements appearing in sect. I. Let us consider an homogeneous 
Pfaffian form in three variables (the generalization to n > 3 variables is straightforward): 

u = P dx + Q dy + Rdz, (Al) 

where P, Q, R are homogeneous functions of degree k (non necessarily k integer) in x, y, z: 

{xd x +yd y +zd z )P = kP- (A2) 
{xd x +yd y +zd z )Q = kQ; (A3) 
{xd x +yd v +zd z )R = k R. (A4) 

The Liouville vector field is Y = (xd x + yd y + zd z ). The integrability condition is 

P(d z Q-d y R) + Q (d x R-d z P) + R(d y P-d x Q) = 0. (A5) 

Then an integrating factor for uj is given by 

f i = Px + Qy + Rz, (A6) 

provided that /i ^ 0. Let us check that 

w e =- (A7) 
A* 

is closed. We have to check the equality of the mixed derivatives, that is, we have to verify that 

d z (-) = d x (-) (A8) 



-)="•(- 
?)"*(f 



dy [ -) . (A10) 



We verify explicitly only (A8). We have 



d z [-) = - d z P - 4 (x d z P + y d z Q + z d z R + R) (All) 
= \{Qyd z P + Rzd z P-Pyd z Q-Pzd z R + PR). (A12) 



Analogously, we get 



d x ( - ) = \ {P x d x R + Q y d x R - R x 8 X P - R y d x Q - P R). (A13) 



We have to show that 



We obtain 



{Qyd z P + Rzd z P~Pyd z Q-P zd z R + P R) 

-(P xd x R + Qyd x R-Rxd x P-Ryd x Q-PR) = 0. (A14) 



y Q {8 X R - 8 Z P) - R {x 8 X P + z 8 Z P) - Ry d y Q 

+P (x 8 X R + z d z R) + P y 8 Z Q = 0. (A15) 



The homogeneity of P, R allows to modify the second and the fourth term above 

yQ(d x R-d z P)+RydyP-kPR-Ryd v Q 
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-Pyd y R + kPR + Pyd z Q = 0. (A16) 

This can rewritten as follows 

y [P (8 Z Q - d y R) + Q (d x R- d z P) + R (d y P - d x Q)\ = 0. (A17) 

Because of the integrability condition (A5), the term appearing in the square parentheses is zero, and the proof is 
completed. This proof is an extension to the case of three variables of the proof appearing in Rcf. [23], p. 19. 
It is also interesting to note that, an analogous generalization of Rcf. [23], p. 18-20, allows to show that, if the coefficients 
P/n, Q/n, R/n of the 1-form uu e appearing in (A7) are homogeneous of degree a / — 1, so that uo e is an exact 1-form 
of degree a + 1 , then the solutions of the equation uo e = are 

g = const. (A18) 

where 

g =-}—(x-+y9. +z ?L\ (A19) 



a + 1 \ /j, fi \i 
is an homogeneous function of degree a + 1 which satisfies 

u> e = dg. (A20) 
In fact, one easily finds that the potential associated with u e is given by 

g = — |— - i Y (w e ) (A21) 
a + 1 

[proof: (a + 1) dg = d iy(u> e ) = — iy du> e + Ly u e = Ly u e = (a + 1) u e }. If, instead, a = — 1, then the solutions of 
uj e = have to be found by quadratures [23]. 



APPENDIX B: FURTHER NOTES ON THE INTEGRAL OF co/fi 

Let us assume that 

n 

cj {k+1) = Y. P * dxl ( B1 ) 

i=l 

is an homogeneous integrablc Pfaffian form of degree k + 1, with Pi homogeneous of degree k for alii = 1 , . . . , n. The 
Liouville operator is 

n a 

y-^w- ,B2 » 

i=l 

If 

n 

H = xi p i (B3) 

i=l 

is the corresponding integrating factor, then 

dW=°^ = d ^ r (B4) 
/i W 

where W > is defined by W = log(M^). Moreover, we know that W has to be an homogeneous function of degree 
q, because dLyW = LydW = 0. We then find 

dw = J2 ( d * w ) dxl = — J2 p * dx ^ ( B5 ) 
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thus one has 



As a consequence, one obtains 



From (B3) one finds 



(O t W) = ^P t Vi = l,...,n. (B6) 



Y W = J2 x ' ( d i W ) ( B7 ) 
»=i 

■ W 

= £ ^ - p i (B8) 

i=l ^ 

= 5 VK. (B9) 



— P 4 = IT (BIO) 



thus the only allowed value for q is q = 1. If g = 0, one would get IT = everywhere, which is absurd. Then IT 
is homogeneous of degree one (extensive function). One then has W(k+i) = dW, where g( fc ) is homogeneous of 
degree k. 

We show that, if G is an homogeneous function of degree q such that 

u (k+i) = 9(k+i- q ) dG, (Bll) 
where g(k+i-q) is homogeneous of degree fc + 1 — g, then necessarily 

G = C W\ (B12) 

where C = const. In fact, one has 

W(fe+i) = 5(fe+i- 9 ) rf G = 5(fc) rf IT ( B1 3) 

and 



Moreover, the homogeneity of G(W) implies 



DG = qG= — W, (B15) 



that is, 



dG G 

^ = *w' (B16) 

whose solution is (B12). 

Notice that q — 1 is a consequence of the construction of the integrating factor [i. In fact, one could find another 
integrating factor p, = fi/q. The corresponding potential would be then G — W q . 

APPENDIX C: IMPLEMENTATION OF ASSUMPTION H8) 

Let us assume that assumption h7) holds and use the variable B = U — b(V, X 1 , . . . ,X n ). We then have V = 
(0, oo ) x Kb- We further assume that the third law (53) holds and that lim^^oo S(U, V, X 1 , . . . , X n ) = +oo for all 
fixed (V, X 1 , . . . ,X n ) £ /Q,. The latter hypothesis is reasonable [notice that it does not hold for systems which allow 
negative temperatures]. Then, for each positive constant c £ IR+, one finds for all fixed (V, X 1 , . . . , X n ) £ JCb 
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lim (S(B,V,X 1 ,...,X n )-c) = -c<0 (CI) 

lim (S(B,V,X 1 ,...,X n )-c) = +oo. (C2) 

Then, a variant of the implicit function theorem [31] shows that there exists a unique function 

B C (V,X 1 ,...,X") : K b -» (0,oo) such that 5(B C (V, A" 1 , . . . , X n ), V, X 1 , . . . , X n ) - c = 0. 

This function is continuous (actually C 2 ) and it defines the unique leaf corresponding to S' _1 (c). 

[Sketch of proof: Let us put Y = V, X 1 , . . . , X n in the following. Then Y £ JCb- S is continuous and monotonically 
strictly increasing in B, hence also the function cr c (B, Y) = S(B, Y)—c has these properties. In particular, a c (B, Y) is 
continuous as a function of B for any fixed Y £ JCb- Permanence of the sign implies that there exist two real numbers 
61 < 62 such that <t c (6i, Y) < < ct c (&2, 50 for all Y £ /Cfc. Then, from the intermediate value theorem it follows that 
there exists a value b £ (i>i, 62) such that a c (b, Y) = for any fixed Y £ JCb- Monotonicity ensures that b is unique. The 
function B C (Y) is then defined as the map B C (Y) = b for each fixed Y £ JC b . This function is also continuous. In fact, 
the set Z(a c ) = {(B ,Y)\a c (B ,Y) — 0} is a closed set, being a c a continuous function. Given a sequence {Y n } C JCb 
such that lim^oo Y n = Y £ JC b , one has (B c (Y n ), Y n ) £ Z(a c ) and, moreover, (lim^oo B c (Y n ), lim^oo Y n ) £ Z(a c ), 
being Z(a c ) closed, lim^oo B c (Y n ) = B C (Y ) then follows.] 
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